The Lepidopterological Society of Japan 


IRER 7y2 to Ga, 38 (4): 269 - 286, 1987 


Species Diversity in the Island Populations of Papilio bianor 
Complex in Japan Considered from Gene Flow and Genetic Distance 


Hirohiko KOMINAMI, Osamu YAMAGUCHI* and Kaoru SUMIYOSHI 


Department of Science Education, Department of Biology* and Department of Applied Physiology, 
Hyogo University of Teacher Education, 
Yashiro-cho, Kato-gun, Hyogo, 673-14 Japan 


Abstract A total of 163 adult males of P. bianor complex were collected in 1986 
at three localities on the mainland of Japan, on the main island of Okinawa, and on 
Ishigaki Island. The butterflies collected from these localities are classified tax- 
onomically into separate subspecies and make up a P. bianor complex. Fluid extract 
of the individual male was analyzed with gel-electrophoresis at six isozyme loci. The 
gene flow and genetic distances between these populations were examined using the 
resultant allelic frequencies. Specific alleles were fixed independently, both in the 
alcohol dehydrogenase locus and in the acid phosphatase locus between the three 
populations. There was no indication of geneflow between the populations. The 
remaining four isozyme loci depicted similar patterns of divergence, but with less 
prominence. Genetic distances, calculated by the formula of NEI (1972), were 0.896 
between populations on the mainland of Japan and on the main island of Okinawa, 
0.737 between populations on the mainland of Japan and on Ishigaki Island, and 1.268 
between populations on the main island of Okinawa and on Ishigaki Island, 
respectively. The magnitudes of these distances strongly suggest that the divergence 
observed in the island populations of P. bianor complex is significantly greater than 
that of interpopulationals within a single species. Allopatric speciation is in progress 
with the aid of random fixation of genes due to the limited number of individuals in 
island populations of P. bianor complex. 


Introduction 


Geographical isolation is an initial stage of speciation. Japanese islands located in 
the middle part of the northern Pacific archipelago from Kamchatka to Formosa seem 
to provide good experimental materials for evolution. There are more than one 
thousand islands, some of which provide suitable habitats to butterflies. Geographical 
isolation is therefore expected to have a significant effect so that some types of 
speciation may be observed. The swallow-tail butterfly, Papilio bianor, is distributed 
widely in the Far East, and it is found throughout the Japanese archipelago. This 
butterfly has a beautifully colored pattern on its wings, which attracts collectors. 
Geographical variations in color pattern and color tone have frequently been observed 
and described (cf. KAWAZOE and WAKABAYASHI, 1976). Numerical taxonomy has 
been applied to distinguish some local varieties (e.g. IZUMI et al., 1984). At present, 
they are divided into seven subspecies, based mainly on morphological differences 
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(FUJIOKA, 1975). Experiments in hybrid inviability and hybrid sterility between 
subspecies have been tried, but the outcome has not yet been sufficient to confirm this 
point because of difficulties in the breeding of larvae. However, data have gradually 
been accumulated (e.g. AE, 1971 and 1985; HAMA, 1976; WAKI, 1981). Recently, 
KAWAZOE and WAKABAYASHI (1976) and ICHINOSE and HORIUCHI (1985) described 
isolates on the Okinawa Islands and on Ishigaki Island as new distinct species derived 
from P. bianor. Taken together, these experimental results suggest that speciation of 
P. bianor complex is proceeding in the northern Pacific archipelago. In this study, an 
experiment was designed to detect genetical differentiation between some island 
populations of this species complex, and to estimate how they are genetically separat- 
ed from each other. Isozyme genes were chosen in this experiment because they are 
considered to be effectively neutral to natural selection and therefore suitable for use 
in the genetic distance analysis (cf. a review of KIMURA, 1982). 


Materials and Methods 


Materials: Adult males of Papilio bianor CRAMER, 1777, were collected at the 
following three localities. (1) Fifty-one males were collected at Akazai, Haga-cho, 
Hyogo prefecture in late May, 1986. These samples represent stocks of subspecies of 
P. bianor dehaani. (2) Fifty-two males were collected at Hanechi, Nago, on the main 
island of Okinawa in early April, 1986. The samples correspond to stocks of subspecies 
of P. bianor ryukyuensis. However, a recent description by KAWAZOE and WAKA- 
BAYASHI (1976) has classified this subspecies into the new sibling species of P. 
okinawensis. Thus, in this article the samples from the main island of Okinawa were 
tentatively named P. okinawensis. (3) Sixty males were collected at a vicinity of Mt. 
Bannadake on Ishigaki Island in late September, 1986. These samples correspond to 
stocks of subspecies of P. bianor okinawensis (=P. bianor junia JORDAN) (KAWAZOE 
and WAKABAYASHI, 1976). 

As a control experiment, Papilio maackii tutanus FENTON, 1881, a sibling species 
to P. bianor, and two unrelated species, Luehdorfia japonica LEECH, 1889, and L. 
puzilot inexpecta SHELJUZHKO, 1913, were employed. Fifty-three adult males of P. 
maackii were collected at Akazai, Haga-cho, Hyogo prefecture in late May, 1986. 
Nineteen males of L. japonica were collected at Nishiwaki, Hyogo prefecture in mid 
April, 1986. Twenty-nine males of L. puziloi were sampled at Fujimicho, Nagano 
prefecture in mid May. These butterflies were immediately put into a deep freezer 
after collection and kept at —80°C. Before use, samples were dissected. The isolated 
testes and arrayed tissues were homogenized with a small amount of distilled water in 
a ceramic grinder. After a short period of centrifugation at 15000 rpm, the supernatant 
was utilized as sample solution. 

Electrophoresis : Isozyme assays were carried out using 5% polyacrylamide gel 
or 1% agarose gel electrophoresis. The concentration and pH condition of gel buffer 
and electrode buffer is shown in Table 1. Running time was 2.0 hours with a current 
of 50 mA. Assayed enzymes were malate dehydrogenase (Mdh), alcohol dehydrogenase 
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Table 1. Buffers and components of stain solution for six enzymes used 
for gene analysis of P. bianor complex. 




















Enzyme Mdh Adh Tox Acph Est-a@ Est-8 
Gel and 90mM EBT buffer 90mM EBT buffer 90mM EBT buffer 13mM Phosphate 13mM Phosphate 13mM Phosphate 
electrophoresis (pH 8.7) (pH 8.7) (pH 8.7) buffer buffer buffer 
buffer (pH 6.8) {pH 6.8) (pH 6.8) 








Stain solution |50mM Tris-HCI 50mM Tris-HCI 50mM Tris-HCI Distilled water Distilled water Distilled water 
buffer (pH 8.5) buffer (pH 8.5) buffer (pH 8.5) 


5mM DL-Malic 2.5% Isopropanol 7mM a@a-Naphthy!] 2% eg-Naphthyl 8% 8-Naphthyl 
acid acid acetate acetate 
phosphate in acetone in acetone 
0.4mM NAD 0.4mM NAD 0.4mM NAD 0.005% Acetic 1mM Naphthanil lmM Naphthanil 
0.3mM NBT 0.3mM NBT 0.3mM NBT acid diazo blue B diazo blue B 
0.4mM PMS 0.4mM PMS 0.4mM PMS 1mM MgCl, 
20mM Sodium 


acetate 
1mM Naphthanil 
diazo blue B 








(Adh), tetrazolium oxidase (Tox), acid phosphatase (Acph), esterase-e (Est-a) and 
esterase-/ (Est-@). The stain method followed STONE et al. (1968) with slight modifi- 
cation, which is also shown in Table 1. The first three enzymes are classified as 
substrate-specific enzymes, and the latter three as substrate-nonspecific ones (KO- 
JIMA, GILLESPIE and TOBARI, 1972). 


Experimental Results 


1. Genetic variability at isozyme loci 

Frequencies of the electrophoretic patterns detected in each population of P. 
bianor are shown in Table 2. The allelic frequencies at each isozyme locus are 
described below. 

Malate dehydrogenase (Mdh): There existed three alleles of Mdh distinguishable 
by agarose gel electrophoresis, the fastest allele toward the cathode (named Mdh”), the 
middle allele (Mdh") and the slowest allele (Mdh‘), in the Akazai population of P. 
bianor. Their frequencies were 0.235, 0.549, and 0.216, respectively. There were three 
alleles in the Ishigaki population of P. bianor, Mdh『, Mdh” and Mdh?, whose fre- 
quencies were 0.267, 0.450, and 0.283, respectively. The mobility of Mdh” was interme- 
diate between Mdh" and Mdh“, and that of Mdh? was intermediate between Mdh™ and 
Mdh*. There were also three segregated alleles, Mdh", Mdh” and Mdh?, at the 
frequencies of 0.288, 0.442 and 0.269 in the Nago population of P. bianor (=P. 
okinawensis). Nomenclature of the alleles is in alphabetical order; ie. the higher 
letter is of faster mobility. The above three populations shared the common allele 
Mdh" at almost the same frequencies. In the control experiment, three alleles Mdhr, 
Mdh® and Mdh” were observed in the Akazai population of P. maackii. Two alterna- 
tive alleles of Mdh? and Mdh? were detected in the Nishiwaki population of L. 
japonica. Monomorphism fixed with Mdh” was observed in the Fujimicho population 
of L. puziloi. These allelic frequencies are summarized in Table 3. An example of the 
zymogram of Mdh is shown in Figure 1. The locus of Mdh was fairly high polymor- 
phic, however a heterozygotic zymograph had never been observed in these materials. 





NII-Electronic Library Service 


The Lepidopterological Society of Japan 


272 Hirohiko KOMINAMI, Osamu YAMAGUCHI and Kaoru SUMIYOSHI 


Table 2. Phenotypic frequencies and their deviation from Hardy-Weiberg equibrium at 6 
isozyme loci in three natural populations of Papilio bianor complex and a 
population of its sibling species of P. maacki and two populations of its unrelated 
species of Luehdorfia japonica and L. puziloi. 


Species P. bianor P. bianor P. okinawensis |P. maackii L. japonica L. puziloi 
Locality AKAZAI ISHIGAKI NAGO AKAZAI NISHIWAKI FUJIMICHO 
No. ind. 51 60 52 53 19 29 
























































Locus SC Exp. ae Obs. Exp. ye OF Exp. ots Exp. nye Oe Exp. P nen O Obs, Exp. 
F 12 120 F 16 160IF 15 150|F 29 290IR 8 80T 29 29.0 
Mdh M 28 28.0/H 27 270IH 23 230|G 23 23.0/T 11 110 
S 11 110|P 17 170|P 14 140IN 1 10 
FF 44 42400 60 600/GG 50 481lEE 8 133100 19 190/00 29 29.0 
FS 5 46 00 2 01lEF 1 00 
SS 0 01 GO 0 38IHH 6 113 
OO 2 01 EH 37 245 
Adh FO 0 36 oO 1 00 
OS 0 02 FO 0 00 
HO 0 10 
EO 0 10 
EF 0 10 
FH 0. 09 
Dev. from | v?ar21 = 0.35 Xf a=2 =11.92** 
H.W. eg. 
SS 51 510|SS 49 495/SS 52 520| FF 6 96ITT 14 127IVV 29 290 
Tox UU 0 05 SS 14 175IVV 2 06 
SU 11 100 FS 33 259ITV 3 57 
Dev. from xar: = 0.03 ¥ate2 =3.99 Vaen=043 | 
H.W. eq. 
FF 10 74IKK 8 48/QQ 52 520ILL 52 510IEE 19 190/DD 29 29.0 
FM 5 1LIÍKP 8 108 00 1 00 
Acph MM 6 41/PP 9 60 LO 0 20 


FS 14 13.0|KR 10 13.6 
MS 12 9.7| PR 12 15.2 














SS 4 5.7|RR 13 96 
Dev. from | x7ar-4= 4.87 Katey = 7.04 
H.W. eq. 
Est-« FF 51 51.0|FF 60 60.0|SS 52 52.0|SS 53 53.0| TT 19 19.0) TT 29 29.0 
FM 5 45/LL 14 18.7| MM 52 52.0| MM 51 51.0|UU 19 19.0|UU 29 29.0 
MM 40 40.6|LR 39 29.6 KM 1 09 
MS 5 4.5|RR 7 11.7 TM 1 0.9 
KM 1 0.9 KK 0 0.0 
Est-8 FF 0 Ol TT 0 0.0 
SS 0 0.1 KT 0 0.0 
KK 0 0.0 
FS 0 0.2 
FK 0 0.1 
KS 0 0.1 
Dev. from | yas-1 =0.04 Xatq2 = 6.06* 
H.W. eq. 


























*Significant at 5% level. 
**Significant at 1% level. 


A possible explanation seems to be that the Mdh locus is sex-linked in Papilionidae 
butterflies (see Discussion). 

Alcohol dehydrogenase (Adh): A small amount of variability was detected in this 
locus through agarose gel electrophoresis. There were three alleles segregated in the 
Akazai population of P. bianor. The main allele was Adh" at a frequency of 0.912. 
The remaining two alleles, Adh? and Adh?, were at the low frequencies of 0.049 and 
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Fig. 1. An example of the zymogram of malate dehydrogenase (Mdh) in the subspecies 
of Papilio bianor complex and its related species of P. maackii. Isozymes of Mdh 
are stained in a deep blue color. White spots are the isozymes of tetrazolium 
oxidase (Tox). P. okinawensis is synonymous with P. b. ryukyuensis collected at 
Nago, Okinawa. P. bianor from Ishigaki is synonymous with P. b. okinawensis 
collected at Ishigakijima, Okinawa. P. bianor from Akazai is synonymous with 
P. b. dehaanii collected at Akazai, Hyogo. (+) and (一 ) stand for the directions 
toward anode and cathode, respectively. 


0.039, respectively. The latter indicated no enzyme activity and was named null allele 
(Adh°). In adult males of the Ishigaki population of P. bianor, no individuals showed 
any Adh activity at all, and were thus scored as monomorphic fixed with Adho. There 
were two alleles Adh* and Adh? at the frequencies of 0.962 and 0.038, respectively, in 
the Nago population of P. okinawensis. Therefore, there was no indication of gene 
flow except the negligible null alleles between these three populations of P. bianor 
complex one with another. In its sibling species, P. maackii, a total of four alleles, 
Adh*, Adh", Adh” and Adh9, were segregated at appreciable frequencies. In both of its 
unrelated species of L. japonica and L. puzilot, no individual males showed any Adh 
activity. Thus, they were scored as monomorphic fixed with Adh?. An example of the 
Adh zymogram is shown in Figure 2. A heterozygote has two separable bands as in 
seen in the figure. Thus, the enzyme seems to be a monomer in protein structure. 
Allelic frequencies are also shown in Table 3. 

Tetrazolium oxidase (Tox): This isozyme analysis was carried out basically on the 
same gel plate as for Mdh or Adh. The locus was rather monomorphic with the 
commom allele Tox? in every population of P. bianor complex. Populations of both 
Akazai and Nago were fixed with this allele. The Ishigaki population contains an 
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Fig. 2. An example of the zymogram of alcohol dehydrogenase (Adh) in the subspecies 
of Papilio bianor complex, in its related species of P. maackii, and in non-related 
species of Luehdorfia japonica and L. puziloi. Isozymes of Adh are stained in a 
deep blue color. For other explanation, see the legend of Figure 1. 
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Fig. 3. An example of the zymogram of tetrazolium oxidase (Tox) in the subspecies of 
Papilio bianor complex. Isozymes of Tox are in white spots. For other explana- 


tion, see the legends of Figures land 2. 
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Table 3. Allelic frequencies and average heterozygosities at 6 isozyme loci in three natural 
populations of Papilio bianor complex, in a population of its sibling species of P. 
maackii and in two populations of its unrelated species of Luehdorfia japonica and 


























L. puzilot. 
Sane P. bianor P. bianor P. okinawensis P. maackii L. japonica L. puziloi 
No- or wenomes AKAZAI ISHIGAKI NAGO AKAZAI NISHIWAKI FUJIMICHO 
6 102 120 104 106 38 58 
examined 
Locus Allele 
Mdh" 0.235 0.267 0.288 0.547 0.0 0.0 
Mdh® 0.0 0.0 0.0 0.434 0.0 0.0 
IMdh『 0.0 0.450 0.442 0.0 0.0 0.0 
Mdh™ _0.549 0.0 0.0 0.0 0.0 0.0 
Mdh Mdh” 0.0 0.0 0.0 0.019 0.0 0.0 
Mdh? 0.0 0.283 0.269 0.0 0.0 0.0 
Mdh? 0.0 0.0 0.0 0.0 0.421 0.0 
Mdhs 0.216 0.0 0.0 0.0 0.0 0.0 
Mdh" 0.0 0.0 0.0 0.0 0.579 1.000 
Adhë 0.0 0.0 0.0 0.500 0.0 0.0 
Adh” 0.912 0.0 0.0 0.019 0.0 0.0 
Adh Adh 0.0 0.0 0.962 0.0 0.0 0.0 
Adh" 0.0 0.0 0.0 0.462 0.0 0.0 
Adh? 0.049 0.0 0.0 0.0 0.0 0.0 
Adh? 0.039 1.000 0.038 0.019 1.000 1.000 
Tox" 0.0 0.0 0.0 0.425 0.0 0.0 
Tox* 1.000 0.908 1.000 0.575 0.0 0.0 
Tox Tox” 0.0 0.0 0.0 0.0 0.816 1.000 
Tox" 0.0 0.092 0.0 0.0 0.0 0.0 
Tox” 0.0 0.0 0.0 0.0 0.184 0.0 
Acph? 0.0 0.0 0.0 0.0 0.0 1.000 
Acph® 0.0 0.0 0.0 0.0 1,000 0.0 
Acphf 0.382 0.0 0.0 0.0 0.0 0.0 
Acph* 0.0 0.283 0.0 0.0 0.0 0.0 
Acphr 0.0 0.0 0.0 0.981 0.0 0.0 
Acph Acph™ 0.284 0.0 0.0 0.0 0.0 0.0 
Acph? 0.0 0.317 0.0 0.0 0.0 0.0 
Acph® 0.0 0.0 1.000 0.0 0.0 0.0 
Acph® 0.0 0.400 0.0 0.0 0.0 0.0 
Acph* 0.333 0.0 0.0 0.0 0.0 0.0 
Acph9 0.0 0.0 0.0 0.019 0.0 0.0 
Est-g『 1.000 1.000 0.0 0.0 0.0 0.0 
Est-@ Est-e* 0.0 0.0 1.000 1.000 0.0 0.0 
Est-a? 0.0 0.0 0.0 0.0 1.000 1.000 
Est- BF 0.049 0.0 0.0 0.0 0.0 0.0 
Est- 8* 0.010 0.0 0.0 0.009 0.0 0.0 
Est- と 0.0 0.558 0.0 0.0 0.0 0.0 
Est-@ Est-8™ 0.892 0.0 1.000 0.981 0.0 0.0 
Est- 8° 0.0 0.442 0.0 0.0 0.0 0.0 
Est- 8° 0.049 0.0 0.0 0.0 0.0 0.0 
Est- 87 0.0 0.0 0.0 0.009 0.0 0.0 
Est- 8” 0.0 0.0 0.0 0.0 1.000 1.000 
Average heterozygosity 0.270 0.328 0.120 0.268 0.131 0.000 
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additional allele Tox” at the low frequency of 0.092. The Akazai population of P. 
maackii also had the common allele Tox’ and another additional allele Tox’. There 
are two alleles, Tox’ and Tox’, in the population of L. japonica. The Fujimicho 
population of L. puziloi had no variability and was fixed with Tox’. An example of 
the zymogram is shown in Figure 3. As shown in the figure, a heterozygote has only 
two distinguishable bands and no additional hybrid bands. Thus, the enzyme is a 
monomer in protein structure. Allelic frequencies are shown in Table 3. 

Acid phosphatase (Acph): Enzyme assay was done by polyacrylamide gel electro- 
phoresis. High variability and high allelic divergence was observed in this locus. 
Every allele had a mobility toward the anode. In the Akazai population of P. bianor, 
three alleles AcphY, Acph™ and Acph* were segregated at almost the same frequecies 
of 0.382, 0.284 and 0.333, respectively. Similarly, in the Ishigaki population of the 
species, there were also three alleles, Acph*, Acph? and Acph®, segregating at the 
approximately same frerquencies of 0.283, 0.317 and 0.400, but each of them had 
different mobility from those of the Akazai population. On the other hand, in the Nago 
population of P. okinawensis, there was no variability in Acph and it was fixed with 
the allele Acph® which had its own mobility. Each allele was endemic to its own 
population. Thus, this may be useful as a good diagnostic characteristic which enables 
us to distinguish separate species one from another. The P. maackii population 
contained another two alleles Acphr and Acph®. Those of L. japonica and L. puziloi 
were fixed with another two alleles, Acph® and Acph?, respectively. An example of the 
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Fig. 4. An example of the zymogram of acid phosphatase (Acph) in the subspecies of 
Papilio bianor complex. Isozymes of Acph are stained in a brown color. Other 
description is same as in Figures 1 and 2. 
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Acph zymogram is represented in Figure 4. The heterozygote carried an additional 
hybrid band to its original two bands. Therefore, the enzyme must be a dimer in 
subunit structure. The situation is easily recognized in a heterozygote between alleles 
AcphF and Acph’ in Figure 4. Allelic frequencies are summarized in Table 3. 

Esterase-a (Est-a): The isozyme patterns were analyzed by agarose gel 
electrophoresis. At the enzyme locus, intrapopulational or intraspecific variation was 
not observed at all. Every population was uniform and fixed with a single specific 
allele. The Akazai and Ishigaki populations of P. bianor were fixed with the common 
allele Est-e『. The Nago population of P. okinawensis and the Akazai population of P. 
maackii were fixed with another common allele Est-e*. Both of the Luehdorfia 
populations were fixed with the other common allele Est-a7. An example of the Est-a 
zymogram is shown in Figure 5. In the figure, esterase-e is one of substrate- 
nonspecific enzymes so that another esterase is also stained on the same gel plate. The 
number of subunits of the enzyme protein was unable to be determined because of the 
non-existence of the allelic heterozygote. 

Esterase-B (Est-@): The isozyme patterns were analyzed by agarose gel 
electrophoresis. A total of four alleles, Est-@*, Est-@*, Est-@™ and Est-@*, were 
segregated in the Akazai population of P. bianor. Their frequencies were 0.049, 0.010, 
0.892 and 0.049, respectively, in this order. In the Ishigaki population of this species, 
two alleles Est-@! and Est-/F were segregated at almost even frequencies of 0.558 and 
0.442. The polymorphism did not overlap with that of the Akazai population. On the 


ESTERASE A 


Genotype. 
【 


Inds No. 





Fig. 5. An example of the zymogram of esterase-a (Est-a) in the subspecies of Papilio 
bianor complex. Isozymes of Est-g are stained in a deep brown color. D. 
melanogaster means the control isozymes of a fruit fly. 
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Fig. 6. An example of the zymogram of esterase-@ (Est-8) in the subspecies of Papilio 
bianor complex. Isozymes of Est-/ are stained in a red color. Other description 
is same as in Figures 1 and 2. 


other hand, there was no Est-/ polymorphism and fixation with Est-@™ in the Nago 
population of P. okinawensis. This allele was also the major one in the Akazai 
population of P. bianor. This was also the case in the P. maackii population, in which 
other two alleles, Est-@* and Est-7, were segregated at low frequencies. In both of 
the Luehdorfia populations, the common allele Est-8V was shared and fixed. An 
example of the zymogram is shown in Figure 6. The allelic heterozygote shows only 
both parental bands and indicates a structural monomer of the protein. These allelic 
frequencies are summarized in Table 3. 

All the above zymographs at 6 isozyme loci are schematically shown in Figure 7. 
Relative mobilities of the isozymes are represented either in cathode direction or in 
anode direction from the origin. 

2. Average heterozygosity at isozyme loci 

Heterozygosity at a single locus (H) is defined by 1- Xlx;’, where x: is the gene 
frequency of the i-th allele. The average heterozygosity (H) is obtained by an 
expansion of the formula into multiple loci. 


a1 Ta 2x”) 


r 


H= 


In this formula, r stands for the number of loci. The quantity is a good indicator of 
the magnitude of genetic variability, and is proportional to the effective size of the 
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Genotype 


MdhF 
Hdh9 
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Mdhs 
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AdhE/AdhE 
AdhF/AdhF 
Adhe/AdhS 
AdhH/AdhH 
AdhS/AdhS 
Adh9/Adh9 
AdhF /AdhSs 
Adh&®/Adh# 


Toxf/ToxFf 
ToxS/Toxs 
ToxT/ToxT 
ToxY/Tox" 
ToxvV/ToxVv 
Tox$/ToxU 
ToxF/Tox$ 
Tox「/ToxVv 


Acph?/Acph? 
AcphE/AcphE 
AcphF /AcphF 
Acph* /Acphk* 
Acpht/Acpht 
Acph"/Acph" 
AcphP/AcphP 
Acph*/Acph®& 
Acph® /Acph® 
AcphS/Acph$ 
Acph0/Acph0 
AcphFf /Acph" 
AcphF /Acphs 
Acph"/Acphs 
AcphK/AcphP 
AcphK/AcphK 
AcphP/AcphR 


Est- cv F/Est-a 
Est-a s/Est- a 
Est-a T/Est-a 


Est-B F/Est- B 
Est-B */Est- 8B 
Est-B '/Est- B 
Est-B "/Est- B 
Est-B ®/Est- B 
Est-8 S/Est-B 
Est-B ™/Est- 8B 
Est-8 Y/Est- 8 
Est-B F/Est-B 
Est- B "/Est- B 
Est-B */Est- B 
Est-B ‘/Est- B 
Est-B '/Est- B 


Fig. 7. Summary of schematic representation of the isozymes at 6 loci examined. 
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population. The average heterozygosity of the Akazai population of P. bianor was 
0.270, that of the Ishigaki population of this species was 0.328, and that of the P. 
okinawensis population was 0.120. These values are not so small such as in a marginal 
population. In the population of its related species of P. maackii, the average hetero- 
zygosity was 0.268. For unrelated species, the average heterozygosity was 0.131 in the 
Nishiwaki population of L. japonica, and no heterozygosity was observed in L. puziloi 
population. This homozygosity seems to be due to the marginal location of the 
Fujimicho, Nagano, population for the distribution range of L. puziloi. Comparison of 
these H values is shown in Table 3. 
3. Genetic distances between the subspecies of P. bianor complex 

Using the above difference in gene frequencies, we can evaluate magnitude of 
genetic differentiation between the subspecies of P. bianor complex. Genetic similar- 
ity (I) and genetic distance (D) are calculated as follows (NEI, 1972). 


ke 
J =X (yu)? 


お <i T 


Jw= > 3 Gury) 
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P.bianor 
AKAZAL 
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Fig. 8. The genetic distances between the subspecies of Papilio bianor complex, its 
related species of P. maackii, and its unrelated species of Luehdorfia japonica and 
L. puzilot. P. bianor from Akazai represents the subspecies of P. b. dehaanii. P., 
bianor from Ishigaki is synonymous with the subspecies of P. b. okinawensis. P. 
okinawensis is synonymous with the subspecies of P. b. ryukyuensis. 
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I= Jxy 
Fy 
D テ ー log。 I 


In the formulae, xi; and yi; mean frequencies of i-th allele of j-th locus in population 
x and population y, respectively. r stands for number of loci examined. 

The Genetic distance between the Akazai and Ishigaki populations of P. bianor 
was 0.737. The distance between the Akazai population of P. bianor and the Nago 
population of P. okinawensis was 0.896. The distance between the Ishigaki population 
of P. bianor and the Nago population of P. okinawensis was 1.268. All of these values 
are significantly greater than the described values of subspecies which range from 
0.028 to 0.201, and are on the same level as that of described full species which range 
from 0.18 to 2.54 with an estimated mean value of 1.0 in insects. This mean value holds 
good also in all described species of animals, plants, fungi and bacteria (a review of 
NEI, 1975). Estimation of genetic distance between another separate species of 
Papilionidae butterflies was also tried in the control experiment. The distance 
between the Akazai populations of P. bianor and P. maackii was 1.009, and the value 
between L. japonica and L. puziloi was 0.241. These are in accordance with the above 
values, and are within the range for full species. The estimated genetic distances 
between the three subspecies of P. bianor complex in question seem to be genetically 
differentiated in the order of species. The genetic distances between all possible 
combinations of the subspecies and the species are illustrated in Figure 8. The genetic 
distance between the genera of Papilio and Luehdorfia was estimated to be 4.229. 


Discussions 


1. The locus of malate dehydrogenase is possibly sex-linked in Papilionidae butter- 
flies 

High genetic variability was observed at the Mdh locus in Papilio species 
examined in the present experiment. The average heterozygosity of the locus was 
0.482. Thus, if we assume that the locus is on one of the autosomes, about a half of 
the number of examined individuals would be expected to show a heterozygotic 
zymo-graph. A total of 264 males were examined, however, none of them showed this 
type of zymograph but a hemizygotic one. A simple explanation is that the locus is on 
one of sex chromosomes. If so, then we can further expect a heterozygotic zymograph 
to be in the female. A total 36 females were examined in the same way, but none of 
them showed heterozygosity, contrary to an expectation. The female also showed 
hemizygotic zymographs as seen in Figure 9. These results taken together mean that 
in both sexes, the enzyme Mdh is produced from one of the genes on alternative 
homogous chromosomes. This type of gene expression is seen only in sex-linked loci. 
On the other hand, the karyotype of butterflies has been thuoght to be heterokryotypic 
in females (ZW typed) based on the low frequency of females among intra- or inter- 
specific hybridization (e.g. AE, 1971). Thus, a possible explanation is that the Mdh 
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Fig. 9. The zymogram of malate dehydrogenase (Mdh) in the female butterflies of two 
species of Papilio. Isozymes of Mdh are stained in a deep blue color. There is no 
heterozygous zymograph formed by heterozygosis. All of them show a hem- 
izygotic pattern. 


locus is located on the Z chromosome in Papilio butterflies, and a mechanism of dosage 
compensation is involved in its expression. Further experiment is necessary to confirm 
this point. 

2. Species diversity of P. bianor complex in island populations 

Many investigators have continuously noticed that there is an extensive variation 
in the morphological characteristics in local populations of P. bianor. Color tone is 
brighter in nothern populations and is darker in southern populations. Some indices of 
wing traits are variable from population to population (IZUMI et al., 1984). The shape 
of male genitalia is highly variable and is unique to some local isolates (KAWAZOE and 
WAKABAYASHI, 1976; ICHINOSE and HORIUCHI, 1985). Interpopulational hybridiza- 
tion has been carried out mainly by AE (1971 and 1985), HAMA (1976) and WAKI (1981), 
disclosing an unbalanced sex ratio in F, progeny and female sterility in some cases. 
From these observation, P. bianor in Japan is subdivided into five subspecies, Papilio 
(Achillides) bianor dehaanii, P. b. amamiensis, P. b. tokaraensis, P. b. ryukyuensis, and 
P. b. okinawensis (FUJIOKA, 1975). However, gene flow between the subspecies in 
nature is still unknown. 

In this experiment, this possibility was directly examined through an analysis of 
isozyme genes. An enzyme is a direct product of a gene. Furthermore, a limited 
number of individuals in a population have a greater effect than the preferential 
pressure of natural selection on isozyme molecules, and make them selectively nearly 
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neutral. Isozyme analysis is therefore suitable for this kind of experiment. At both of 
the Adh and Acph loci, unique alleles were segregated independently between the 
populations of P. b. dehaanii, P. b. ryukyuensis and P. b. okinawensis as seen in Table 3. 
There was no indication of gene flow between the populations. On the other hand, 
migration of butterflies between islands has been observed frequentlly (e.g. AZUMA 
and MINATO, 1983). Sterility of F, hybrid females from a cross between P. b. ryukyu- 
ensis and P. b. dehaanii was reported by HAMA (1976). This is also the case in a cross 
between P. b. okinawensis and P. b. dehaanii (WAKI, 1981). Reproductive isolation 
between the subspecies seems to be perfect in nature. 

Genetic divergence between the subspecies was also surveyed quantitatively in 
terms of the genetic distance, and the results are given in Figure 8. The genetic 
distances in any combinations between the subspecies were greater than 0.7. The 
magnitude was of the same order as that between P. b. dehaanii and P. maackii 
tutanus, or was more than that between L. japonica and L. puziloi. This magnitude is 


L. japonica 0.241 L.puziloi 


NISHIWAKI FUJIMICHO 
i 
4.229 
i 
i 
P.bianor | 1.009 P.maackili 
AKAZAT AKAZAI 
0.737 | ^ 0.896 
P.bianor FF . P.okinawensis 
ISHIGAKI Än ~~ NAGO 
1.268 


Fig. 10. The genetic distances between the subspecies of Papilio bianor complex 
(triangle), its related species of P. maackii, and its unrelated species of Luehdorfia 
japonica and L. puziloi. P. bianor from Akazai represents the subspecies of P. b. 
dehaanii. P. bianor from Ishigaki is synonymous with the subspecies of P. b. 
okinawen--sis. 

P. okinawensis is synonymous with the subspecies of P. b. ryukyuensis. 
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also significantly greater than the mean value of genetic distances between subspecies 
examined so far, in the order of 0.1, and does not fall within the range of those of 
conspecific local populations, in the order of 0.01 (a review of NEI, 1975). Recent work 
by MATSUOKA et al. (1984) calculated the genetic distance between a pair of mor- 
phologically very similar species of Neope goschkevitschii and N. niponica to be 0.753. 
The summary of the present results is shown in Figure 10. 

Adaptation of a subspecies to its own local habitat is possibly estimated through 
the amount of average heterozygosity (H). The value of H of P. b. ryukyuensis was 
0.120, and that of P. b. okinawensis was 0.328, as seen in Table 2. Both of these did not 
so deviate too much from that of P. b. dehaanii, 0.270. Thus, each subspecies popula- 
tion has comparable genetic variability, and it is further suggested that it had a long 
time, to build up to the amount of present genetic variabililty after the foundation of 
the ancestral population. 

The above observations are strongly suggestive of speciation in progress in Papilio 
bianor complex, especially in island populations. Some of present subspecies of P. 
bianor may be reclassified into new independent species. Recently, KAWAZOE and 
WAKABAYASHI (1976) described subspecies of P. b. ryukyuensis as a new species of P. 
okinawensis based on morphological differences. Similarly, ICHINOSE and HORIUCHI 
(1985) described subspecies of P. b. okinawensis as a new separate species based on 
genetical analysis. The present experimental results support these author’s conclu- 
sions and give additional experimental evidence. Similar speciation is thought to be in 
progress also in other butterflies in Japanese southwest island populations. 
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摘 要 


ヵ カラ ス ア ゲハ と その 近 縁 種 間 の ァ アイ ッ ザイ ム 人 遺伝 子 頻 度 よ り 解 析 し た 分 化 程 度 
(小南 裕彦 ・ 山 口 修 ・ 住 吉 薫 ) 


カラ ス ア ゲ へ (P. bianor complex) は , 中 国 大 陸 ・ 朝 鮮 ・ 台 湾 ・ 沖 縄 諸島 ・ 日 本 本 土 に 広く 分 布 し , その 
は ね の 模様 ・ 色 合い も 各 分 布地 で 微妙 に 変化 し て いる . この 傾向 は 特に , 南西 諸島 の よう な 海 で 隔 た っ た 隔 
離 小 集団 で は 強い が , 形態 学 的 に は 現在 7 亜 種 に 分 類 さ れる こと が 多い . 川副 ・ 若 林 (1975) で は , この うち 
沖縄 諸島 産 を P. okinawensis と し , 一 ノ 瀬 ・ 堀 内 (1985) は , 八重 山 諸島 産 の も の を 別種 と し て いる . 

し か し , これ ら 亜 種 間 で の お 互い の 交配 実験 も 多数 試み られ て は いる が , F。 ま で の 飼育 が 困難 で も あり , 
種 (species) レ ベル まで の 分 化 が 進ん で いる の か , ある い は さら に 多数 の 別種 の レベ ル ま で 分 化し て いる の 
か , 実験 的 な 確証 は 得 ら れん て いな い . 

本 研究 で は , これ ら の 関係 の 一 端 を 調べ る た め , 1986 年 に カラ ス ア ゲハ の 本 州 産 , 沖縄 本 島 産 , 石垣 島 産 
の 雄 を 各 51, 52, 60 RSL, 

腹部 を 解剖 し , 精巣 お よび 消化 器官 に 含ま れる 酵素 の うち , リン ゴ 酸 脱水 素 酵 素 (Mdh), アル ュー ル 脱 水 
素 酵素 (Adh), テト ラ ゾ リ ウム 酸化 酵素 (Tox), MERRIER (Acph), ヶ - エ ステ ル 化 酵素 (Est-a), 8-5 
ステ ル 化 酵素 (Est-8) の 6 遺伝 子 座 に つき , その 対立 遺伝 子 頻度 を 電気 泳動 法 で 調査 し た . 同一 の 移動 ペ 
ター ン お よび 特異 的 移動 パコ ター ン を 持つ も の の 割合 より , NEI(1972) の 式 を 用 いて , 互い の 遺伝 的 距離 を 推 
定 し た . 遺伝 的 距離 は 2 集団 間 で の すべ て の 同一 移動 度 を 持つ アイ ソ ザ イム 遺伝 子 の 頻度 の 相関 の 逆数 と 
し て 表 さ れる . 

本 州 産 と 沖縄 本 島 産 の ぁゃ の で の 遺伝 的 距離 は , 0.896 で あり , 本 州 産 と 石垣 島 産 , 沖縄 本 島 産 と 石垣 島 産 
の それ は , 0.737 と 1.268 で あっ た . これ ら の 値 は , 対照 区 に と っ た カラ ス ア ゲハ と ミヤ マカ ラス アゲ ハ ( ア 
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maackii) D 1.009, BLO, ギフ チョ ウ (L. japonica) と ヒメ ギフ チョ ウ (L. puziloi) の 0.241 と 同 程度 の 
オー ダー 万分 化し て いる と と を 示し て いる . これ ら の 値 は , NEI の 総論 (1975) の 他 の 動植物 で 得 ら れ て い 
る 種 間 の 値 ; 0.18 - 2.54 の 範囲 内 に 入っ て いる . 

同一 種 内 の 他方 品種 間 で は , 0.01 の オー ダー, BRHICH, 0.1 の オー ダー, 属 間 で は , 1.10 以上 の オー 
ダー か ら み て も 種 の レベ ル に 最も 近い . 

また , Adh お よび Acph で は , 3 産地 で 共に 共通 の 対立 遺伝 子 は 全く 存在 し て お ら ず , 野外 で は 遺伝 子 の 
交流 が 無い こと を 示す と と も に , 固有 の 対立 遺伝 子 に 固定 し て いる こと を 示し て いる . 

これ ら の こと は , 南西 諸島 の カラ ス ア ゲハ は , 遺伝 的 に は 別種 と し て も 良い 程度 に 分 化 を 起こ し て いる こ 
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